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Large scale traffic networks are an indispensable part of contemporary human mobility and inter- 
national trade. Networks of airport travel or cargo ships movements are invaluable for the under- 
standing of human mobility patterns[l], epidemic spreading[2], global trade[3] and spread of invasive 
species [4]. Universal features of such networks are necessary ingredients of their description and can 
point to important mechanisms of their formation. Different studies[5] point to the universal char- 
acter of some of the exponents measured in such networks. Here we show that exponents which 
relate i) the strength of nodes to their degree and ii) weights of links to degrees of nodes that they 
connect have a geometric origin. We present a simple robust model which exhibits the observed 
power laws and relates exponents to the dimensionality of 2D space in which traffic networks are 
embedded. The model is studied both analytically and in simulations and the conditions which 
result with previously reported exponents are clearly explained. We show that the relation between 
weight strength and degree is s{k) ~ k^^'^, the relation between distance strength and degree is 
s''(fc) ~ k^^'^ and the relation between weight of link and degrees of linked nodes is vUij {ki kj ) 
on the plane 2D surface. We further analyse the influence of spherical geometry, relevant for the 
whole planet, on exact values of these exponents. Our model predicts that these exponents should 
be found in future studies of port networks and impose constraints on more refined models of port 
networks. 
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Study 


/3 




a 


Barrat et al.[6] World Airport Network (WAN) 


1.5 ±0.1 




0.5 ±0.1 


Barrat et al.[7] North Amorican Airport Network 


1.7 


1.4 




Bagler [8] Airport Network of India 


1.43 ±0.06 






Zhang et al.[9] 


1.58 






Passenger Airport Network of China 








Zhang et al.[9] 


2.2 






Cargo Airport Network of China 








O'Danleyman et al.[W] (WAN) 


1.33 






Kaluza et al. [11] Network of Global Cargo Ship 
Movements (NGCSM) 


1.46 ±0.1 






Hu and Zhu[12] World Wide Maritime Trans- 
portation Network 


1.3 






O'Danleyman et al. [10] (NGCSM) 


1.33 







TABLE I. Measured exponents in different studies. Except for the exponent /? of Cargo Airport Network of China, all the 
reported exponents fit well to the predictions of our model. Li and Cai[13] also suggested the value of a = 1/2 for the Chinese 
Airport Network as an ansatz for their data, although without clear fitting procedure. Most of the studies did not report error 
intervals for the measured exponents but, based on those that did, we can assume errors of around 0.1. 



An intense empirical research of traffic networks of aeroplanes[l, 6-10, 13, 14] and cargo ships[ll, 12, 15] points 
to the fact that there are exponents of power laws universally present in these types of networks. In general, nodes 
in such traffic networks represent ports. Weighted links of these networks represent fluxes of some units - people 
or goods - between these ports [6]. In general these fluxes are directed, but it is also common to represent overall 
fluxes between ports as undirected, which is the approach that we follow in this study. The weighted hnk Wjj found 
in empirical networks therefore represents some average number of transported objects back and forth, during some 
time interval between ports i and j. The strength Si of the port i is calculated as a sum of fluxes Wij of units through 
the port i.e. Si = J2j ^ij- The distance strength of the port sf is a sum of distances to the neighbouring ports i.e. 

= J2jen(i) ^iJ' where n{i) represents the set of ports connected to the port i. Different studies [1, 6, 7, 11-15] 
observed interesting relations between aforementioned properties and topological variables. Namely, the weighted 
strength of the port i depends on its number of neighbours ki as Si ^ fcf [1, 6, 11-15]. The distance strength of the 
port i scales with its number of neighbours ki as sf ~ fef"* [7]. Finally the flux between two ports of degrees ki and 
kj scales as Wij ^ [kikj)'^ [6, 13]. Values of these exponents in different studies are presented in Table I. Exponents 
in airport networks have been much more studied than exponents of cargo ship networks and we do not know for any 
attempt of measuring a in cargo ship networks. Nevertheless, agreement between findings in cargo ship networks and 
predictions of our model lead us to believe that our model could capture the behaviour of cargo ship networks as well. 

Although a number of models coupling topological and traffic properties of such spatial networks has been proposed, 
none of these models managed to explain all three observed exponents. In[16] authors proposed a model based on 
spatial preferential attachment [17]. They found a relation between exponents, which explains only one or the other 
observed exponent, but not all of them. Further, in[7, 18] authors studied the model of weight driven attachment which 
was able to explain some of the observed properties of transport networks, but failed to reproduce exponents which 
relate degree properties with weight properties. Their model yields the value of /? = 1, and the authors proposed that 
the measured exponents are a consequence of some non-linear process lying behind the attachment model they have 
used. The most complete analysis of coupling between topology and weights is laid out in[19], in which the author 
proposed two models - a weighted attachment model with the addition of weights and a fitness based model[20]. These 
models can reproduce any possible value of /3, including /3 = 3/2, and a number of other network properties, but are 
sensitive to tunable model parameters, and exponents Z?'' and a were not reproduced. 

We observe that values of measured exponents point to the possibility that they are rational numbers. Traditional 
statistical mechanics [21] suggest that rational exponents can often be attributed to the dimensionality of space. If the 
values of measured exponents are really of geometric origin, then a fairly general model should be able to reproduce 
these exponents. Such a model should produce relatively stable exponents with respect to other model parameters 
and possible additional refinements of the model. 

We propose a simple and robust mathematical model of port networks. This model is based on three assumptions. 
First, we assume that the probability for an individual unit to travel from one port to the other is determined by 
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some hidden variables assigned to all ports. These variables depend on properties of the port such as the size of 
population, its economic power, number of international companies, etc., but could also depend on the traffic network 
funcionality to include transient ports etc. The predictions of the model do not depend on the precise definition 
of these variables or the precise mechanism of their assignement to ports. These hidden variables are modelled as 
the fitness of every port - a random number x, which is attached to every port and is drawn from some probability 
distribution p{x) [20] . We will later show that exponents are very robust to the choice of probability distribution and 
to possible spatial correlations between different fitness variables. The second assumption is that the probability of 
intended travel between two ports is proportional to the product of their fitness variables p{i — )■ j) ~ XiXj. This 
assumption is related to the usual model of random mixing in networks[22]. The fitness variables represent some 
units which travel back and forth to interact with some other units. The probability that some unit from port i 
will travel to interact with some unit from port j is, in its simplest instance, proportional to the amount of units in 
port j. This assumption means that the expected flux Wij of units from port i to port j and from port j to port i 
is proportional to the product of fitness variables The third assumption is related to the economic and 

geometric properties of traffic networks. We assume that the link between ports i and j will exist only if the expected 
volume of traffic will cover the costs of distant travel. The costs of travel Cij from port i to port j are modelled as a 
linearly increasing function of their distance r^ . This assumption is related to the fact that costs of travel are, among 
other things, determined by the amount of fuel one has to spend travelling from one place to another and by the 
costs of crews which are also proportional to the distance via the number of travels a crew can manage in the unit of 
time. We also present the data on the travel fares between different U.S. airports and their mutual distances [23]. The 
relationship between distance and fares has Pearson's correlation r = 0.64 (Supplementary Figure 1), which justifies 
our assumption of linearity in the model. Notice that this data state real commercial fares which are not equal to the 
lower bound of transport profitability, which is relevant for this model. Nevertheless, the linear relationship between 
fares and distance is still clearly present. We expect that this relationship would be even more pronounced in the 
data on pure costs of travel between different ports. 

The proposed model is studied both analytically and in simulations. Let us first consider an idealized case in which 
ports are randomly spread on an infinite 2-dimensional plane. We assume that the distribution of ports is relatively 
homogeneous and model their spatial distribution as a spatial Poisson process [25]. To every existing port i a random 
number X{i) drawn from some probability distribution p[x) is assigned. Then a relation between the expected degree 
{k{X)) of the port and its fitness variable X can be written as 

{k{X)) = a J dA e{f{xX)-c{r))p{x)dx. (1) 

Here ct is a surface density of ports, dA is a differential element of the surface, is a Heaviside step function, / 
is function of earnings generated from the flux of travelling units, and c is a cost function of travel which grows 
with distance. We assume that these functions are monotonously increasing. In the rest of this paper we assume 
f{u) = c(it) = M as the simplest linear relationship. In Supplementary Information 1 we give a short analysis of 
some other possibilities for functions / and c which lead to the same exponents. In it, we also show that the relation 
needs to be satisfied on the average only to produce the same exponents as the exact version of the model. 
Similarly to equation (1) the expected strength and expected distance strength of the port can be related to its 
fitness X as 



{s{X)) J dAj xXQ{xX -r) p{x)dx. (2) 
(X)) = a J dAj^ rQ{xX- r) p{x)dx. (3) 



is' 



The integrations in these equations arc easily performed to obtain {k{X)) — ■KaM2{p)X , {s{X)) ~ T:aM^{p)X^ and 
{s'^{X)) — ^'kgM^{p)X^ . Here M2^z{p) represent the second and the third moment of the distribution p{x). From 
these relations immediately follows that {s{k)) ~ k^/"^, {s'^{k)) ^ k^l"^ and {w{ki,kj)) ^ {{ki){kj)y/^. For simplicity 
we refer to these exponents as infinite exponents in the remaining text. One can also understand this relationship 
qualitatively. Port i will be connected with some port j up to certain radius rf^ given by fitness values of these ports 
X{i) and X{j). Some port i' which has the double fitness of the port i, X{i') = 2X{i) will connect with port j on 
double distance r^,j = 2r|j-. Since we also assume the homogeneous distribution of ports, the number of ports to 
which port i' is connected grows with the surface i.e. k ^ ^ X^ . Weights are assumed to be proportional to the 
fitness X and to the number of connections k which yields s ~ Xr^ ~ X^ . The average distance of the port is also 
proportional to the distance and to the typical number of connections s'^ ^ rk ^ X^ . Having this in mind, it is clear 
that the relationship between the cut off distance r'^ and the fitness X has to be linear in order to reproduce measured 
exponents. 
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Although distribution p{x) can, in principle, be freely chosen on the infinite 2D plane, it is constrained for the 
applications by the observed degree distributions Po{ko) extracted from the data. If we assume that ko{x) is well 
approximated by {k{x)) and since the relationship between expected degree {k{x)) and the fitness variable x is a 
monotonic function k ~ J-(x) ^ x^ , then the relation p{x) ^ \^^^^\po{T{x)) follows. In the case of scale free 
distributions Po{k) ~ k"'' , the relation yields p{x) ~ x'^~'^'' . Relevant studies[6, 8, 12, 15] have found power law 
distributions of degrees in port networks and in validations of our model we investigate power law distributions of 
fitnesses for exponents 76 (1,5]. Another constraint on the distribution is the requirement that every port should be 
reachable from every other port. In other words, we demand the existence of a giant component in the network[24]. 
Yet another constraint is that the network should be sparse enough as the real port networks are. This condition 
can be controlled by monitoring the average degree of the network. All these factors have to be taken in account to 
choose model parameters, since not every distribution of fitnesses will satisfy all constraints for all densities or for all 
functions / and c. These constraints were accounted for by considering the minimal and maximal values of fitness 
variables. 

In Figure 1 we show dependences of weighted properties on topological properties for finite 2D space depending on 
the exponent 7 of the power law distribution p{x). For a broad range of power law exponents simulated exponents 
a, P and P'^ are close to theoretically predicted values, as can be seen in Figure 1. The increase of the observed 
exponents as power law is becoming steeper which is related to the process of network sparsification, as shown in 
Figure 2. In the limiting case of a network consisting of disconnected subgraphs it is hard to talk about exponents at 
all. The degree distribution and the aforementioned exponents are clearly not the only properties of port networks 
previously reported. In Supplementary Information 2, for the reasons of completeness, we also present an analytical 
calculation for the clustering coefficient in infinite 2D plane. For thelargc k we observe behavior close to C{k) ^ k~^, 
a feature previously reported by other researchers [6, 7]. 

In the Supplementary Information 3 we present an analytical treatment of the model for the finite 2D plane and 
for the sphere of radius R. The case of the sphere is particularly important since it represents the real geometry in 
which port networks operate. The expected degree {k{X)), the strength {s{X)) and the distance strength {s'^{X)) of 
the port with fitness X on a sphere can be computed from relations 

27T TT 00 

{k{X)) = aR'^ J dip J sinOde J p{x)dx, (4) 
Re/x 



{s{X)) ^uR"- [ dip I sin0d0 J Xxp{x)dx, (5) 

R.e/x 



and 

27T TT OO 

(s'^iX)) ^ crR^ J dip J sinOde J eRp{x)dx. (6) 

Re/x 

The change of the geometry influences the dependences of these variable on the fitness X, and we find deviations 
from the infinite exponents. Nevertheless, the average distances between ports reported in the previous studies, point 
to the fact that the majority of the traffic takes place on the locally almost fiat elements of a sphere. The average 
distance {d) of travels in the case of the North American airport network is around 1000 km[7], which means that the 
relative difference between sin{{d) / R) and {d)/R is only around 0.4%. We believe that this is the main reason why 
the monitored exponents of port networks really resemble theoretical values calculated for the infinite 2D plane. 

In real traffic networks, ports are not homogeneously distributed. In the previous text we dealt with density of 
ports which arc scattered over the space as Possionian random variables [2 5]. To check if the heterogeneities which 
necessarily exist among the real ports infiucncc the behaviour of exponents, we gathered the data on geographical 
locations in the U.S. [26]. In Figure 2 we present simulations of our model on the geographical locations of American 
airports. 

In Figure 3, we present three different sets of data points taken from the database [26] and construct a network among 
them. The exponents do not differ much among these different spatial port distributions, as long as the distribution 
of the fitness variable fulfils the aforementioned constraints. The robustness of this model to the heterogeneity of 
spatial port distribution is clear. It is also clear that many ports in the world are used just as transit airports for a 
significant amount of traffic and the fitness variables used in the model do not capture such behaviour. On the other 



5 



hand, such usage of airports also influences the realistic possible values of fitnesses and we expect that in reality they 
have a certain spatial hierarchical distribution. 

In conclusion, in this paper we presented a very robust model which can be altered in many conceivable ways, 
which generally reproduce observed exponents, which relate weighted properties to topological properties. Even 
further it can be shown that in the case of general D-dimensional flat space, the exponents would be a = \/D and 
jS ^ 13'^ ^ [D + 1)/D (see Supplementary Information 4), which could (apart from intellectual curiosity) possibly be 
interesting for the case of interstellar traffic with D=3. Similar exponents of geometric origin have been found in the 
case of allometric scaling in biological systems [27], and at present it is not clear if the presented model is related 
to these phenomena. Furthermore, we believe that some other modes of transportation like, for instance, intercity 
buses, could also be well described with this model, but at present we are not aware of any additional research in 
that direction. Since a convincing simple relationship between fitness variables proposed in this paper and measurable 
socio-economic quantities is still lacking, further research in this direction is clearly needed. 



I. METHODS 

All simulations of the model were performed in FortranQO. The locations of the ports were drawn randomly from 
uniform distribution defined on the circle of area 10 in arbitrary units in the case of a finite 2D plane. We used 
different truncated power-law distribution for fitness variables. The parameters of the distribution were: minimal 
value of fitness Xmin = 0.1, maximal value of fitness Xmax = 1-6 and exponent 1 < 7 < 5. For every set of points 
or different parameters we have performed 100 different simulations with stochastic variables drawn from power-law 
distributions. In the case of the sphere with radius R = 1 we have used exactly the same parameters. The locations of 
ports on the real maps were acquired from geographical data. We have used geolocations of of airports, hospitals and 
capes in the U.S. [26] in order to encapsulate as many of real spatial nonpoissonian heterogeneous point distributions 
as possible. For the Geolocations we have used Xmin = 0.01 or Xmin = 0.02 and Xmax = 0.65. We have also used 20 
realizations, due to a large number of data points. The exponents were calculated by the least square fitting procedure 
of the logarithms of related variables. Reported errors are standard deviations of fitted exponents in the simulation 
sample. 
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SUPPLEMENTARY INFORMATION 1 



The cost function is generally defined as ,f{xiX2) — c(ri2). The expected degree of node with fitness X is 

{k{X)) j dA j e {f{xX) - c{r)) p{x)dx 

poo .c-^fixX)) 

= 2'Ka I p{x)dx I rdr 

jQ Jo 

poo 

= 7rcr / [c-^{f{xX))fp{x)dx. 



Similarly the expected strength of the node with fitness X is 

POO 

{s{X))=na / [c-\f{xX))f{xX)p{x)dx. 
Jo 

Equations (7) and (8) again give the same exponent behaviour if we choose f = c 

pOO 

{k{X))^na [r\f{xX))fp{x)dx 
Jo 

= 7raX^M2{p). 

{S{X))=7T<7 / [r\f(xX))]'(xX)p{x)dx 

Jo 

= iraX^Msip). 



(7) 



(8) 



(9) 
(10) 



In the case of different functions we can again easily reconstruct proposed exponents. For instance let f{x) = fx, and 
c(x) — c\x + C2. Then equations for the s and k are: 



k{X) 
s{X) 



pXH'hip) Jc^XAhip) , ci 



pX^Nhip) ^fc2XHh{p) , 4XAh{p) 



2fc2XAh{p)-ci 



^ 1 and 



2fc2X^M2ip)-c:iXMi{p) 



(11) 

^ 1 . In practice 



Proposed exponents will be reconstructed as long as 

it will workj allways when fX{X) ^ C2. 

Let us further consider the case in which flow Wij is not necessarily exactly equal to the product XiXj. Realistically, 
one can expect that some ports are better connected than expected by chance, while some others are not so well 
connected. For example, airports of two neighbouring states which do not have diplomatic relationships because of 
political reasons can have much lower traffic than one would expect from fits of their overall traffic. This additional 
preference of the links can be modelled with a new random variable rj taken from some distribution p{r]). Let us assume 
that the average value {rj) = I. Traffic flows can be written as Wij = rjijXiXj. If r] and x are independent random num- 
bers, the equations for degree, strength and strength distance for infinite 2D plane are {k{X)) = TTaM2{p)M2{p)X'^ , 
{s{X)) = t:(jM^{p)M^{p)X'^ and {s{XY) = 2-K(jM-i{p)M^{p)X^ /i. Obviously, infinite exponents are preserved even 
if we allow some fluctuations around expected values of flows. 



SUPPLEMENTRAY INFORMATION 2 

The expected degree dependent clustering (C(fc)) coefficient of some port on an inffiiite 2D plane can be calculated 



as: 

oo oo 2-7T 2tt oo 

^2 



cr^ / ridri J r2dr2 J dtpi J dip2 J p{xi)dxi J p{x2)dx2<d{Xxi - ri)Q{Xx2 - r2)<d{xiX2 - ri2) 

{C{X)) = ^ ^ °- (12) 

0-2 J ndri J r2dr2 J dfi J dip2 J p{xi)dxi J p{x2)dx2Q{Xxi — ri)Q{Xx2 ~ r2) 





where ri2 = \/rj + — 2rir2 cos{ipi - ip2) 
The denominator D is equal to: 



Xxi oo XX2 



D^Aa'^n'^ J p{xi)dxi J ridpi j p{x2)dx2 J r2dr2 (13) 



(14) 





1 

2^'^2(p) 



and therefore 

D = a^TT'^X^[M2{p)f (15) 

The numerator A is equal to 

oo Xxi oo Xx2 2-iv 2tt 

A^a'^ J p{xi)dxi j ridri J p{x2)dx2 J r2dr2 J d(pi J dip2'c>{xiX2 - ri2). (16) 



We can make a convenient substitution 

^(<<52 - </3l) (17) 
^'-^(<P2+</5l). (18) 

Separating the integral domain in two parts we can write integrals over new angles: 

2TV+ip IT 2Tr — ifi 

d^ [ d^'^-i^:^Q{x,X2~rr2)+ fd^ I d^'^-^^^Q{x,X2 - r,2). (19) 

— 7r —ip 

Using Eq. (19)the numerator can be written as 

oo 

A ^ a"^ j p{xi)dxi j ridri j p{x2)dx2 j r2dr22 j dip{2TT - ip)Q{xiX2 - ri2) (20) 



The random variables are positive semidefinite and thus 



6(^1X2 - Jrl +rl- 2rir2 cos ^) = e(cos ^ - d±Il i^l^y (21) 

2rir2 



We can write 

Xxi 00 XX2 



A = J pix,)dxi J rrdri J p{x2)dx2 J r2dr2 X (^^^^4^^^^^) ' ^^^^ 
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Function x(a;) is defined as 



1 a; < -1 

x(a^) = { ^arccosx — 1 < a; < 1 (23) 
a; > 1 



Finally the clustering coefficient is 

Xxi OO XX2 



{C(X)) = J p{x,)dx, J ndn J P{x2)dx2 J r,dr, x ^ ^J^'"''^' ) (24) 



Behavior of the clustering coefficient with respect to k is presented in supplementary figure 2. 
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SUPPLEMENTARY INFORMATION 3 



Let us consider a point on the pole of the sphere of radius R. The sphere is covered with random points at locations 
9 and 0. The density of these points is given with a — S{r ~ ri)dr/S. The distance of points from the pole 

is given with RO, where R is the radius of the sphere and 9 is the inclination angle. We assume that the point on 
the pole is a port with the variable X and that other points represent different ports whose fitnesses are drawn from 
some unspecified distribution p{x). The expected degree of the port on the pole is then 



{k{X)) =0"^ y J sin9d9 J p{x)dx. 
Re/x 



(25) 



Similarly the strength {s) and distance strength (5^) of the vertex can be calculated from: 

2-K TT 00 



{s{X)) ^<^R^ j j sm9d9 J Xxp{x)dx, 







Re/x 



(26) 



{s'^iX)) ^^R^ j dip j sin9d9 J 9Rp{x)dx. 



RM/X 



The expected degree can be written as two terms equation as 



7r 00 

k{X) = 2TraR^ J sin9d9 J p{x)dx, 



eR/x 

ttR/X 



2ti(jR^ 



= 2TiaR^ 







1 — cos ■ 



Xx 



p{x)dx + 2 / p{x)dx 



ttR/X 



/ I 1 ^ cos — — 1 p[x)dx + / 

Jo \ R J JttR/X 



1 + cos -j^ ) p(x)dx 



(27) 



and similar equations can be written for strength and distance strength also. The second integral can be neglected 
for distributions whose tail falls fast enough. In the case of our simulations {X^ax = 1-6, i? = 1), the second integral 
is exactly since it lower limit is tt/1.6 « 2 while the distribution p{x > 1.6) = 0. After the calculation of space 
integrals we make an expansion of trigonometric functions up to the second order: 



(fc(X)) = iraX^M^ip) 



12R'^ 



XHh{p) = TT<jXHl2{p) (1 + Cfc) 



(28) 



We can now calculate: 



(s(X)) - naXHMp) - ^xHhip) = naX^Mp) (1 + c,) , 



(s'iX)) = ^iraX^M.ip) - ^X'M,{p) = '^TTaX^M,{p) (1 + c,.) 



^.2..XM,(p)-i^X3M,(p) 



(29) 
(30) 

(31) 



(32) 
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Thus 



d{s}_ 
d{k) 



Sttct 
12i?2 



With the assumption that corrections are smah j^^i^j^j^^'p) ^ 1 and -^"^^;f^;^^^^^^ ^ 1 we get: 



d{s) SA'hip)^ Mi{p) 



X - 



d{k) 2M2{p) 4i?2M2(p) 
Solving (28) for X and expanding in series: 



5M5(p) M3(p) 



(i(s)_ ^ 3 Af3(p) ;^^^i/2 



6M4(p) M2(/9). 

M3(p) _^J^M5(p)■ 



48i?2A^25/2(p)(7^C^)3/2 
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Integrating with condition that (s) = when (fc) = wc get: 

Miip) 



Similarly 



Ah{p) 



^Ml"{p) 
M^ip) 



(fc) 



3/2 



-(fc) 



1/2 



24i?2Af|/^(p)(7rr7)3/2 

Af4(p) 



M2(p) Nh{p) 
Msip) ^M,{p) 



{k) 



3/2 



d{k) ^mI''^{p) 24i?2M|/2(p)(7ra)3/2 



M2(p) M4(p) 



\ M^{p) _ ^M5(p) 



(fc) 



5/2 



Ahip) AU[p) 



{k) 



3/2 



(33) 



(34) 



(35) 



(36) 



(37) 



3 ^mI'^{p) 60RHl^J^{p){7ra)3/^ 

Leading order corrections to degree Ck and strength Cg are negligible if i2H^mI{p) ^ ^ ^^^^ 12H^mI{p) ^ ^' ^® 
calculated l.h.s. of this relations for the simulations described in this paper. These contributions were computed for 
the minimal and maximal exponent of p{x) and for mutual couplings between ports with maximal fitness value as 
well as couplings between ports with maximal fitness value and ports with typical fitness value. 

The influence of the spherical geometry in the case of the sphere with homegnous distribution of ports can be 
quantified via typical values of the flows i.e. products of fltnesses: X,nin = 0.1 X„iax ~ 1-6: 

• X^ = Xl^^, 7 = 1 : Cfc = 0.274, c, - 0.328 

. X2 = X^aAX), 7 = 1 : Cfe = 0.0927, = 0.1106 

• ^2 = X2„^^, 7 = 2.6 : Cfc = 0.136, c, 0.2296 

• X2 = X,naAX): 7 - 2.6 : Cfc - 0.0185, c, ^ 0.0314 
United states Airports: X,nin ~ 0.02, Xmax = 0.65: 



.^J^ MP) 



5/2 



(38) 



• X2 = X2„^^, 7 = 1 : Cfe = 7.445 • lO^^, c, = 8.926 • lO^^ 

• X2 = X„,aAX), 7 = 1 : Cfc = 2.073 • 10-^, = 2.485 • lO^^ 

• X2 = X2„„^, 7 = 2.6 : Ck = 3.298 • lO^^, c, = 6.172 • lO^^ 

• X2 = X,naAX), 7 = 2.6 : Ck = 2.380 • 10-^, = 4.454 • lO'^ 

As shown, the influence of spherical geometry is practically non-existent for typical situations and can have some 
effect only on ports with fltnesses drawn from the tail of the distribution. Spherical geometry clearly influences the 
relation between strength s, distance strength s'^ and degree k. In general we can not expect the inflnite exponents 
on spherical geometries. Nevertheless, in cases in which most of the links between ports are below a certain threshold, 
the geometry is effectively almost flat and the predicted exponents are expected to be found. 
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Here we calculate the expected values of degree and strength for given X in the case of finite flat geometry. Now 
the domain of integration is finite and, if we assume the convexity of the domain, we can write the equation for the 
expected degree of the port with fitness X as: 



27r r{ip) oo 

k{X) = J dtp J rdr J p{x)dx, 

r/X 



(39) 



where T{(p) describes the distance of the domain edge to the origin with respect to angle ip. This integral can be 
separated in two parts as: 



k{X) = (T J dp 





Ix^ 



riv)/x 



x^p{x)dx + -r^{(p) j p{x)dx 



27r cjo 



■naX"^ j p{x)dx — '2^'^^ j j 



T{v)/X 



For the strength the equation is very similar: 



Xx 



x^ p{x)dx. 



(40) 



(41) 



2tv 00 



s{X) = ttcjX^ j x^p{x)dx - icrX^ j dtp j 



1 



T(^)/X 

These expressions can be rewritten in a more transparent form as: 



Xx 



x^p{x)dx. 



(42) 



and 



00 

'2 / ^2 , 



k{X) = TTCrX^y X^p{x)dx 




00 

-2 I 2 . 



TTcrX J X p{x)dx 




00 

2iT J 

1 f r(^}/x 

2i r^" — 



x'^p{x)dx 



J x^p{x)dx 




I- 



r(^) 

X 



(43) 



(44) 



00 

-3 f ^3, 



s{X) = TiaX* j x-'p{x)dx 




1 - 



00 

-3 13, 



TToX-^ j X p{x)dx 





h 1^ 



x^p{x)dx 



J x^p{x)dx 




1- e 



X 



(45) 



(46) 



If we assume an upper cut-off of the distribution Xmax, an assumption which is very realistic for any imaginable 
port system, then if V(/3 : T(ip) > XXmax both e functions are equal to zero. If this condition is not fulfilled, a 



weaker condition will still preserve infinite limit exponents. Namely, a factor 1 — 



EM 

Xx 



can achieve values from 



for X — ^-^^ up to 1 for large X which reduces the contribution of the second term. Finally, by averaging over the 
angle ip with the assumption of a convex domain the contribution of terms different than zero is reduced by factor 2 
at worst. Finally, if we assume Xmax < where F represents the characteristic linear dimension of the domain, we 
can easily expect the infinite limit exponents recovered too. 



13 

SUPPLEMENTARY INFORMATION 4 

111 the hypothetical case of a D-diniensional infinite plane, our model is integrated over D-dimensional volume with 
the differential element dVu = r^^^C]jdr. From this trivially follows 

/•OO /"OO 

{k(X))^aCD r^-^dr Q {xX - r) p{x)dx, (47) 
Jq Jq 

/•OO poo 

{s{X))^aCD r^'^dr xXQ {xX ^ r) p{x)dx, (48) 
Jo Ja 

/•OO pQC 

{s'^{X))^aCD r^-^dr rQ {xX - r) p{x)dx. (49) 
Jo Jo 

(50) 

One can see that the exponents (3 = fi''- — {D + 1) / D and a — l/D 
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FIG. 1. Simulations for a finite 2D plane are represented in these figures. The upper left figure represents the a dependence on 
the power law exponent 7, the upper right figure represents the /3 dependence on the power law exponent 7, and the middle left 
figure represent the dependence on the power law exponent 7. It is clear that infinite exponents are well reproduced in finite 
fiat space as well. The middle right figure represents the relation between strength and distance strength in the simulated data. 
The red line is the theoretically predicted result of s = 2/3s'*. The lower left figure and the lower right figure represent the 
strength s and degree strength s"^ dependence on degree k. The data points are simulations and the red line has an exponent 
3/2 and is drawn for comparison. Strength s, distance strength s'* and degree fc as a function of fitness X are in the insets. 
The data in insets are exponentially binned to improve visibility and lines with exponents 3 and 2 are drawn for comparison. 
The networks for this simulation are realized on a finite fiat space as described in the Methods section. The exponents are well 
reproduced for a broad range of distribution parameters. 
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1.60 




FIG. 2. The upper left figure represents the q dependence on the power law exponent 7, the upper right figure represents the 
j3 dependence on the power law exponent 7, and the middle left figure represent the /3'' dependence on the power law exponent 
7. Locations of points are locations of 18745(26] American airports and fitness variables are randomly drawn from the power 
law distributions described with exponent 7, minimal value of fitness Xmin ~ 0.02 and maximal value of fitness Xmax = 0.65. 
It is clear that infinite exponents are well reproduced in spherical geometry for 7 < 3. The middle right figure represents the 
dependence of /3 on the fraction of airports P{k > 0) which have at least one neighbour. It is clear that the simulated exponents 
differ from infinite ones when the network starts to disconnect due to a big number of airports with small fitness parameters. 
The lower left figure represents the relationship between /3 and the maximal degree of the airport in the network. Again, when 
the maximal degree is too small, networks effectively disconnect and simulated exponents deviate from the infinite case. The 
lower right figure demonstrates that higher density of points reproduces the infinite exponents better. The relationship between 
/3 and 7 for three cases of the different density of points in the case of a sphere of radius i? = 1 is presented. A'' is the number 
of points on the sphere and the density is cr = N/An. 
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FIG. 3. The upper left figure represents locations of 18745 American airports, the upper right figure represents locations of 
13385 American hospitals and the lower left figure represents locations of 11938 American capes. Colours represent different 
local densities of locations in order to emphasize differences among heterogeneities of this data. The relationship between /3 
and P{k > 0) for these sets of ports is shown in the lower right figure. The minimal value of fitness is Xmin ~ 0.01 and maximal 
value of fitness is Xmax = 0.65. The infinite exponent is again reproduced for all three sets of spatial points, as long as p{x) is 
sufficiently broad i.e. a network is connected. 




FIG. 4. The average fares of travel among 200 largest american airports are placed on the y-axis, and on the x-axis is the 
mutual distance of airports. Data points represent the average fare of travel in the third quarter of 2010 between two airports at 
a given distance. The data was downloaded from [23]. The straight line represents linear regression obtained via least squares. 
Pearson's correlation of the data is 0.64. It is visible that the smallest fares at a given distance also have a linear like behaviour, 
and we can expect that the smallest fares are a better representation of distance costs then the average fares. 
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FIG. 5. An average value of the clustering coefficient C is placed Onon the y-axis, and the square of the fitness is placed 
on the X-axis. The data points (pluses) represent the expected clustering dependent on the square of ^ k. The full line 
represents the power law with exponent —1 and such assymptotic relationship is found in real networks. The minimal and 
maximal value of x are again 0.1 and 1.6, which were used in simulations as well. 



